Thermodynamics of MHD flows with axial symmetry 



N. LeprovoslQ and B. Dubrulle 

DRECAM/SPEC/CEA Saclay, and CNRS (URA2464), F-91190 Gif sur Yvette Cedex, France 

P-H. Chavanis 

Laboratoire de Physique Theorique (UMR 5152), Universite Paul Sabatier, Toulouse, France 

We present strategies based upon extremization principles, in the case of the axisymmetric equa- 
tions of magnetohydrodynamics (MHD). We study the equilibrium shape by using a minimum 
energy principle under the constraints of the MHD axisymmetric equations. We also propose a 
numerical algorithm based on a maximum energy dissipation principle to compute in a consistent 
way the equilibrium states. Then, we develop the statistical mechanics of such flows and recover 
the same equilibrium states giving a justification of the minimum energy principle. We find that 
fluctuations obey a Gaussian shape and we make the link between the conservation of the Casimirs 
on the coarse-grained scale and the process of energy dissipation. 

PACS numbers: 05.70.Ln,05.90.+m,47.10.+g,52.30.-q 



I. INTRODUCTION 

The recent success of two experimental fluid dynamos 0, has renewed the interest in the mechanism of dynamo 
saturation, and, thus, of equilibrium configurations in MHD. At the present time, there is no general theory to tackle 
this problem, besides dimensional theory. For example, in a conducting fluid with typical velocity V , density p , 
Reynolds number Re and magnetic Prandtl number Pm, the typical level of magnetic field reached at saturation is 
necessarily 

B 2 = p oP V 2 f(Re,Pm), (1) 

where / is a priori an arbitrary function of Re and Pm. Many numerical simulations Q lead to / = 1, i.e. equipartition 
between the magnetic and turbulent energy. This is therefore often taken as a working tool in astrophysical or 
geophysical application. However, this result is far from applying to any saturated dynamo. Moreover, it does not 
give any information about possible anisotropy of the saturated field. It would therefore be interesting to build robust 
algorithms to derive the function /. By robust, we mean algorithm which depends on characteristic global quantities 
of the system (like total energy) but not necessarily on small-scale dissipation, or boundary conditions. 

An interesting candidate in this regards is provided by statistical mechanics. In the case of pure fluid mechanics, 
statistical mechanics has mainly been developed within the frame of Eulcr equation for a two-dimensional perfect 
fluid. Onsager [j| first used a Hamiltonian model of point vortices. Within this framework, turbulence is a state of 
negative temperature leading to the coalescence of vortices of same sign |(| ■ Further improvement were provided by 
Miller et al. Q and Robert and Sommeria who independently introduced a discretization of the vorticity in a 
certain number of levels to account for the continuous nature of vorticity. Using the maximum entropy formalism of 
statistical mechanics 0, it is then possible to give the shape of the (meta)-equilibrium solution of Euler's equation as 
well as the fine-grained fluctuations around it 01- This is similar to Lynden-Bell's theory of violent relaxation 
in stellar dynamics (see Chavanis 0] f° r a description of the analogy between 2D vortices and stellar systems). The 
predictive power of the statistical theory is however limited by the existence of an infinite number of constants (the 
Casimirs) which precludes the finding of an universal Gibbs state. In particular, the metaequilibrium state strongly 
depends on the details of the initial condition. In certain occasions, for instance when the flow is forced at small 
scales, it may be more relevant to fix a prior distribution of vorticity fluctuations instead of the Casimirs 0. Then, 
the coarse-grained flow maximizes a "generalized" entropy functional determined by the prior distribution of vorticity 
[l4|. The statistical mechanics of MHD flows has been recently explored by Jordan and Turkington 0] in 2D. In 
contrast with non-magnetized 2D hydrodynamics they obtained a universal Gaussian shape for the fluctuations. This 
comes from the fact that the conserved quantity in the MHD case is an integral quantity of the primitive velocity and 
magnetic fields and thus, in the continuum limit, has vanishing fluctuations. 
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The pure 2D situation however seldom applies to astrophysical or geophysical flows. In this respect, it is interesting 
to develop statistical mechanics of systems closer to natural situations, albeit sufficiently simple so that the already 
well tested recipes of statistical mechanics apply. These requirements are met by flows with axial symmetry. Most 
natural objects are rotating, selecting this peculiar symmetry. Moreover, upon shifting from 2D to axi-symmetric 
flows, one mainly shifts from a translation invariance along one axis, towards a rotation invariance along one axis. 
Apart from important physical consequences which need to be taken into account (for example conservation of angular 
momentum instead of vorticity or momentum, curvature terms), this induces a similarity between the two systems 
which enables a natural adaptation of the 2D case to the axisymmetric case. This is shown in the present paper, 
where we recover the Gaussian shape of the fluctuations and make the link between the conservation of the Casimirs 
on the coarse-grained scale and the process of energy dissipation. 

In the first part of the paper, we study the equilibrium shape by using a minimum energy principle under the 
constraints of the MHD axisymmetric equations. We also propose a numerical algorithm based on a maximum energy 
dissipation principle to compute in a consistent way the equilibrium states. This is similar to the relaxation equation 
proposed by Chavanis 0] in 2D hydrodynamics to construct stable stationary solutions of the Euler equation by 
maximizing the production of a ff-function. Then, we develop the statistical mechanics of such flows and recover 
these equilibrium states, thereby providing a physical justification for the minimum energy principle. 

II. MHD FLOWS WITH AXIAL SYMMETRY 
A. Equations and notations 

Consider the ideal incompressible MHD equations: 

<9 t U + (U-V)U = --VP + (V x B) x B , 
P 

<9 t B + (U-V)B = (B-V)U, (2) 

where U is the fluid velocity, P is the pressure, ^/p/io B is the magnetic field and p is the (constant) fluid density. In 
the axisymmetric case we consider, it is convenient to introduce the poloidal/toroidal decomposition for the fields U 
and B: 

U = U p + Ut =U P + Ue e , (3) 
B = B p + Bt = V x (Ae e )+Be e , 

where A = A p + A eg is the potential vector. This decomposition will be used in our statistical mechanics approach. 

When considering energy methods, we shall introduce alternate fields, built upon the poloidal and toroidal decom- 
position. They are : a u = rU, a b = rA, = oj/r and £j = B/r, where u> is the toroidal part of the vorticity field. 
In these variables, the ideal incompressible MHD equations J5J become, in the axisymmetric approximation, a set of 
four scalar equations: 

d t a b + {tp,a b } = 0, (4) 
d t a u + {ip,a u } = {er 6 ,2y£ b } , 

2 

dtt» + {i/>,tu} = <9 z (^-4 2 )-{a b ,A*a b }, 

where the fields are function of the axial coordinate z and the modified radial coordinate y — r 2 /2 and ip is a stream 
function: U p = V x (ifj/r eg). We have introduced a Poisson Bracket: {f,g} — d y fd z g — d z fd y g. We also defined 
a pseudo Laplacian in the new coordinates: 

d 2 1 d 2 

A * = d^ + 2y"d^- (5) 

Following jlfij . we will make an intensive use of the operators (for more details, see appendix^): curl which gives the 
toroidal part of the curl of any vector and Curl which takes a toroidal field as argument and returns the poloidal 
part of the curl. If j = curlH is the toroidal part of the current and ip = r Cw/ _1 (U P ), the following relations hold: 

e u = -a*v, (6) 

j/r = -A*o-fc. (7) 
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Under the shape J2J| , the ideal axisymmetric MHD equation of motion lead to the immediate identification of a b = rA 
as a conserved quantity associated to axial symmetry. This quantity is only advected by the velocity field and thus 
should play a special role regarding the global conserved quantities, as we now show. 



B. Conservation laws 



1. General case 

The whole set of conservation laws of the axisymmetric ideal MHD equations have been derived by Woltjer [Tr| : 

I = J C{a b )dydz, (8) 

H m = 2 j £ b N(a b )dydz, 

H c = J {F{a b )i u + a u i b F\a b )}dydz, 

L = J a u G(a b ) dydz , 



where C, N, F and G are arbitrary functions. One can check that these integrals are indeed constants of motion by 
using Q and the following boundary conditions: a b = a u = £ u = £ b = on the frontier of the domain. To prove the 
constancy of the third integral, one has to suppose that F(0) = 0. The interpretation of these integrals of motion is 
easier when considering a special case, introduced by Chandrasekhar ]vf\ - 



2. Chandrasekhar model 



The conservation laws take a simpler shape when one considers only linear and quadratic conservation laws, such 
that N(a b ) = F(a b ) = G(a b ) = a b and N(cr b ) = G(a b ) — 1 . The case F(<j b ) = 1 is forbidden by the condition 
that F should vanish at the origin. In that case, the set of conserved quantities can be split in two families. The first 
one is made-up with conserved quantities of the ideal MHD system, irrespectively of the geometry: 



H m = 2 J £ b a b dydz — yA Bdx = 2y AB c?x , (9) 
H c = J{<r b £ u + o- u 6} dydz = J U • B dx. , 

E = \ J je^ - a b A,a b + g + 2y£ b 2 J dydz = ~ J (U 2 + B 2 ) dx , 

where H m is the magnetic helicity, H c is the cross-helicity and E is the total energy. Note that due to the Lorentz 
force, the kinetic helicity is not conserved, unlike in the pure hydrodynamical case. The other family of conserved 
quantities is made of the particular integrals of motion which appear due to axisymmetry: 



C{a b )dydz = J C(rA)dx, (10) 
H' m = 2 J{, b dydz = J ^-dx, 
L = J cT u G{cr b )dydz = J r 2 UBdx, 

L' = J ' <r u dydz = J rU dx . 

Apart from L' the angular momentum, it is difficult to give the other quantities any physical interpretation. The 
class of invariant / are called the Casimirs of the system (if one defines a non canonical bracket for the Hamiltonian 
system, they commute, in the bracket sense, will all other functionals) . The conservation laws found by Woltjer are 
then generalization of these quantities. 
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C. Dynamical stability 

1. General case 

Following we show that the extremization of energy at fixed /, if m , H c and L determines the general form 
of stationary solutions of the MHD equations. We argue that the solutions that minimize the energy are nonlinear ly 
dynamically stable for the inviscid equations. 

To make the minimization, we first note that each integral is equivalent to an infinite set of constraints. Following 
Woltjer, we shall introduce a complete set of functions and label these functions and the corresponding integrals with 
an index n. Then, introducing Lagrange multipliers for each constraint, to first order, the variational problem takes 
the form: 

SE + Yl{a {n) 5I^ + n$8HW + n^SHW + 7 (™)<SL(») j = . (11) 

71=1 ^ ' 

Taking the variations on a b , £&, o u and £„, we find: 

A»<7 6 - -F'(a b )A^ + F"(a b )a u ^ + G'(a b )a u +2N'(a b )^ b + C'(a b ), (12) 
2y& = -2N(a b ) - F'(a b )a u , 

£ = -F'(<7 6 )6-G(a 6 ), 
i/f = -F(a b ), 

where we have set F(a b ) = Yln=o Fn(o~b) an d similar notations for the other functions. This is the general 
solution of the incompressible axisymmetric ideal MHD problem In the general case, it is possible to express 

the three field a u , and in terms of a b . Then the first equation of the above system leads a partial differential 
equation for a b to be solved to find the equilibrium distribution. Note that the extremization of the "free energy" 
J = E + al + [i m Hm + [i-cH-c + yields the same equations as (|12|) . Differences will appear on the second order 
variations (see below). 



2. Chandrasekhar model 



In the Chandrasekhar model, the arbitrary functions arc at most linear functions of o b : N(a b ) = n m <J b + /i' r 
F(a b ) = fi c a b and G(a b ) — 7C7{, + 7'. Thus the stationary profile in the Chandrasekhar model is given by: 



m • 



A*cr h = - /UcA*^ + 7<7 U + 2/z m £b + C'(a b ) , (13) 

2y&> = -2/i m cr fc - 2fJ,' m - fi c C7 u , 
~Hc£,b - jcr b - 7' , 



ijj = -[i c a b . 



From the previous equations, we obtain 

2y(l-A*c)6 = 2(7/icy-Mm)<Tfe + 2Mc7'2/-2/x^, (14) 

(1 - n 2 c )a u = 2(fi c n m - jy)a b + 2^ m - 27'y , 

ip = -n c a b , 



where a b is given by the differential equation: 



(1 - £) 2 ^a b = $(a b ) - [2^ + 2j 2 y]a b - 2^'y . (15) 

y y 

These expressions can be used to prove that these fields are stationary solutions of the axisymmetric MHD equations. 
We now turn to the stability problem. Since the free energy J = E + al + ji m H m + [i c Mc + lL is conserved by the 
ideal dynamics, a minimum of J will be nonlinear ly dynamically stable (at least in the formal sense of Holm et al. 
[l8|). Note that this implication is not trivial because the system under study is dimensionally infinite. We will admit 
that their analysis can be generalized to the axisymmetric case. Since the integrals which appear in the free energy 
are conserved individually, a minimum of energy at fixed other constraints also determines a nonlinearly dynamically 
stable stationary solution of the MHD equations. This second stability criterion is stronger than the first (it includes 
it). We shall not prove these results, nor write the second order variations, here. 
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D. Numerical algorithm 



1. General case 



It is usually difficult to solve directly the system of equations l)14 |l -l)f 5 [l and make sure that they yield a stable 
stationary solution of the MHD equations. Instead, we shall propose a set of relaxation equations which minimize 
the energy while conserving any other integral of motion. This permits to construct solutions of the system (|14|l - (|15|) 
which are energy minima and respect the other constraints. A physical justification of this precedure linked to the 
dissipation of energy will be given in Sec. IIII C II 

Our relaxation equations can be written under the generic form 



9(7 _ , 

where a stands for a b , £&, a u or £„. Using straightforward integration by parts, we then get: 



/ = 



L 

E 



3 crb -[VC'(a b )]dydz, 
2 / |j Cb • V[N(a h )} + 3 ab ■ V[N'(a b )t b ]\dydz, 

J Cu • V[F(a b )} + 3 ab ■ V[F'(a b )^ u + F"{a b )a u ^ b ] 
J CT „ • V[F'(<7&)&] + 3 ib ■ W[F'{a b )a u ]\dydz , 
J ffu • V[G(cr fe )] + J CTb ■ W[G'(<j b )a u ] \dydz , 



(16) 



(17) 



3 iu • W- J CTb ■ V(A*<7 6 ) + J CTu • V^g ] + J., • YC2!/0.)}'1H<I-: ■ 



To construct the optimal currents, we rely on a procedure of maximization of the rate of dissipation of energy E very 
similar to the procedure of maximum entropy production principle (MEPP) of Robert and Sommeria [2lT | in the 2D 
turbulence case. This is equivalent to say that the evolution towards the equilibrium state (|14|) - (|15fl is very rapid. 
We thus try to maximize E given the conservation of /, H m , H c and L. Such maximization can only have solution 
for bounded currents (if not, the fastest evolution is for infinite currents). Therefore, we also impose a bound on 
where, as before, a stands for a b , a u , £ u . 

Writing the variational problem under the form 



+ °° ( ~] 1 



(18) 



and taking variations on 3 ab , 3^ b , J CTu , J^ u , we obtain the optimal currents. Inserting their expressions in the relaxation 
equations (|16f) . we get: 



^ = V-|^ CTb V- -A^ b + C'{a b ,t) + 2£, b N'{cj b ,t) 
+ ZuF'(a b ,t)+a u &F"(a b ,t)+G'(a b ,t)a u 



(19) 



dt 
dcr u 
dt 

dt 



2y^ b + 2N(a b ,t) + F'(a b ,t)a, 



^+^ b F'(a b ,t)+G(a bl t) 



V-S.D 6u W-[i; + F(a b ,t)]\. 



where we have set F(a b , t) — J2n=o l I c n \t)F n (a b ) and similar notations for the other functions. The time evolution of 



the Lagrange multipliers //i (t) etc. are obtained by substituting the optimal currents in the constraints He' 1 ' = etc 



(n) 
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and solving the resulting set of algebraic equations. Using the expression of the optimal currents and the condition 
that I — H m = H c = L = 0, we can show that: 

E -l{k + k + k + lk} iyd *^ m 

provided that the diffusion currents Dg u , D Ub ,-D CTll and D^ b are positive. Thus, the energy decreases until all the 
currents vanish. In that case, we obtain the static equations <|12|) . In addition, this numerical algorithm guarantees 
that only energy minima (not maxima or saddle points) are reached. Note that if we fix the Lagrange multipliers 
instead of the constraints, the foregoing relaxation equations lead to a stationary state which minimizes the free energy 
J. Then, as stated above, the constructed solutions will be nonlinearly dynamical stable solution of the MHD set of 
equations. However, not allowing the Lagrange multiplier to depend on time, we may "miss" some stable solutions of 
the problem. Indeed, we know that minima of the free energy are nonlinearly stable solutions of the problem but we 
do not know if they are the only ones: some solutions can be minima of E at fixed /, H m , H c and L while they are 
not minima of J = E + al + u. m H m + fj, c H c + jL. 



2. Chandrasekhar model 



A j-A*cr 6 + C'(a b , t) + 2n m {t)& + n c {t)i u + j(t)a u } , (21) 
A<^ 2y£ b + 2^ m (t)a b + fi c (t)<J u 



A{V> + fJ, c (t)(Xb}. 



In the Chandrasekhar model (with \x' m = 7' = 0), the previous equations can be simplified. The equilibrium 
solution does not depend on the particular value of the diffusion coefficients (these are only multiplicative factors of 
the optimal currents) and for simplicity, we set = D ab = D au = D^ b = 1. The relaxation equations then reduce 
to: 

d<r b 
dt 

dt 

^ = A{g+ Mc (t)& + 7 (t)o* 
dt 

where the Lagrange multipliers evolve in time so as to conserve the constraints l|17l) . 

These equations are the MHD counterpart of the relaxation equations proposed by Chavanis for 2D hydrody- 
namical flows described by the Euler equation. In this context, a stable stationary solution of the Euler equation 
maximizes a H-function (playing the role of a generalized entropy) at fixed energy and circulation. A justification of 
this procedure, linked to the increase of H-functions on the coarse-grained scale, will be further discussed in Sec. IIVI 
and compared with the MHD case. 

If we set the velocity field to zero (a u — £ u = 0), we get a system of equations linking the poloidal part (a b ) and 
the toroidal part (£&) of the magnetic field. It is fairly easy to see that the coupling between the two quantities is 
proportional to n m the Lagrange multiplier associated to the conservation of magnetic helicity. This is reminiscent 
of the a effect of dynamo theory (see Steenbeck et al. m the "kinematic approximation" where the effect of 

the Lorentz force is removed, the coupling between the toroidal and poloidal part of the magnetic field is given by a 
coefficient proportional to the kinetic helicity of the fluctuating velocity field. Our model is not able to recover this 
fact because, as noticed above, this quantity is not conserved in the full MHD case. However, taking into account the 
retroaction of the magnetic field on the velocity field, Pouquet et al. [2(| were able to write the non-linear a-effect as 
a difference between the kinetic and the magnetic helicity of the fluctuations: a = Hk — H m . Our relaxation equations 
therefore recover the fact that the approach to saturation of the magnetic field is mainly monitored by the magnetic 
helicity. 



III. STATISTICAL MECHANICS OF AXISYMMETRIC FLOWS 



In the previous section, we obtained general equilibrium velocity and magnetic field profiles through minimization 
of the energy under constraints. In the present section, we derive velocity and magnetic field distribution using a 
thermodynamical approach, based upon a statistical mechanics of axisymmetric MHD flows. As we later check, the 
distribution we find are such that their mean fields obey the equilibrium profiles found by energy minimization. For 
simplicity, we focus here on the Chandrasekhar model. 
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A. Definitions and formalism 



Following 0, and 0, we introduce a coarse-graining procedure through the consideration of a length-scale 
under which the details of the fields are irrelevant. The microstates are defined in terms of all the microscopic possible 
fields u(x) and b(x). On this phase space, we define the probability density p(r, u,b) of a given microstate. The 
macrostates are then defined in terms of fields observed on the coarse-grained scale. The mean field (denoted by a 
bar) is determined by the following relations: 

U(x) = J up(r,u,b)dudb, (22) 

B(x) = J b p(r, u, b) dudh . 



We introduce the mixing entropy 



S[p] = - J p(r, u, b) in[p(r, u, b)] drdudb , (23) 



which has the form of Shanon's entropy in information theory j2^| 0. The most probable states are the field U and 
B which maximize the entropy subject to the constraints. The mathematical ground for such a procedure is that an 
overwhelming majority of all the possible microstates with the correct values for the constants of motion will be close 
to this state (see Q for a precise definition of the neighborhood of a macrostate and the proof of this concentration 
property). Note that this approach gives not only the coarse-grained field (U, B) but also the fluctuations around it 
through the distribution p(r, u,b). 

Each conserved quantity has a numerical value which can be calculated given the initial condition, or from the 
detailed knowledge of the fine-grained fields. The integrals calculated with the coarse-grained quantities are not 
necessarily conserved because part of the integral of motion can go into fine-grained fluctuations (as we shall see, 
this is the case for the energy in MHD flows). This induces a distinction between two classes of conserved quantities, 
according to their behavior through coarse-graining. Those which are not affected are called robust, whereas the other 
one are called fragiles. 

B. Constraints 

In this section, it is convenient to come back to the original velocity and magnetic fields. The constraints are the 
coarse-grained values of the conserved quantities ©■ The key-point, as noted by 0, is that the quantity coming 
from a spatial integration of one of the field u or b, is smooth. In our case, it amounts to neglecting the fluctuations 
of A which is spatially integrated from B and write A — A. Thus, the coarse-grained values of the conserved quantity 
are given by: 



1=1 C(rA)dx, (24) 
H„, = 2 / A B dx , 



H c = J u • b p(r, u, b) dxdudb , 
E = ^ J {u 2 + b 2 ) p(r,u,b) dxdudb , 



H' m = 2 I — dx, (25) 



r 

.2 



L = J AUr z dx , 
11 = f Ur dx. 



The constraint / is the Casimir, connected to the conservation of en, along the motions. In the present case, it is a 
robust quantity as it is conserved on the coarse-grained scale. As stated previously, the quantities H m , H c and E are 
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the mean values of the usual quadratic invariants of ideal MHD, namely the magnetic helicity, the cross-helicity and 
the energy. On the contrary, the quantities H' m , L and U are specific to axisymmetric systems. Because these last 
three conservation laws are usually disregarded in classical MHD theory, it is interesting in the sequel to separate the 
study in two cases, according to which the conservation of H' m , L and L' is physically relevant ("rotating case") or is 
not physically relevant ( "classical case" ) . 

C. Gibbs state 

1. Classical case 

The MHD equations develop a mixing process leading to a metacquilibrium state on the coarse-grained scale. It 
is obtained by maximizing the mixing entropy S[p] with respect to the distribution p at fixed I, H m , H c and E (we 
omit the bars in the following). We have: 

SS = - J {I + In p) 5p dxdudb , (26) 

SH C = J u • b Sp dxdudb , 

SE = 1 J(u 2 + b 2 ) 6 p dxdudb. 

The variation of the magnetic helicity and the Casimirs is more tedious because they involve the coarse-grained field 
A. For the magnetic helicity, we have: 



5H m = 2 J (SAB + A 5B)dx . (27) 
Now, using an integration by parts, it is straightforward to show that 

5 A B dx = [ SB P -Apdx. (28) 



Therefore, 

5H m = 2 J (SB P ■ A P + A SB) dx = 2 J A ■ SB dx (29) 
= 2 J A ■ b Sp dxdudb . 
Regarding the variation of the Casimirs, we find: 

SI = J C'{rA)r5Adx = j C'{rA)r Curr 1 B P dx= J curl" 1 [rC'{r A)] ■ 6B P dx, 



(30) 



or 



61 = J curl" 1 [rC'{r A)] ■ b P Sp dxdudb . (31) 
Writing the variational principle in the form 

+oo 

SS - PSE - p rn SH m - p c 5H c - a (n) <y/ (n) = , (32) 



n=l 

we find that 

P r..2 i u2\ r,.. a u .. u l-lr^-t'/ 



1 + lnp = -tt(u 2 + b z ) - 2p m A ■ b p c u ■ b curl" L [rC'(r A)] ■ b P . (33) 
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It is appropriate to write u = U + u' and b = B + b' where the first term denotes the coarse-grained field. Then, the 
equation l|33|) can be rewritten 

1+lnp = ~^(u' 2 + b 12 ) - Mc u' • b' - MmA • B - • B (34) 
- (H + u ') • [PIJ + fi c B] 



(? + b ') ' \P& + 2m™A + nJj + curl" 1 [rC'(r A)} 



Hence the fluctuations are Gaussian: 



P= \ exp|-|(u /2 + b /2 )-At c u'-b'| = iexpji^^^ j , (35) 

where we defined a 6-dimensionnal vector: Xi = (u'i, u' 2 , u' 3 , b' ll b' 2 , b' 3 ). The mean-field is given by: 

0U + n c B = , (36) 
j3B + 2(j, m A + nJJ = , 
0B P + 2fi m A P + n c \J P + curP 1 [rC'{rA)\ = . 

Taking the curl of these relations and using curlBp = j, curlUp = u> and curlAp = B, we recover the equilibrium 
distribution (|13|l with 7 = 7' = p! rn = fi' c = 0. Therefore, in this classical case, the equilibrium profiles are such 
that mean velocity and mean magnetic field are aligned. This is a well known feature of turbulent MHD, which has 
been observed in the solar wind (wh ere v w ±B). It has been linked with a principle of minimum energy at constant 
cross- helicity (see chapter 7.3 of (24J and references therein). This feature is also present in numericals simulation of 
decaying 2D MHD turbulence, where the current and the vorticity are seen to be very much equal |25| . This can 
therefore be seen as the mere outcome of conservation of quadratic integral of motions, and may provide an interesting 
general rule about dynamo saturation in systems where these quadratic constraints are dominant. 

Using the Gaussian shape for the fluctuations, it is quite easy to derive the mean properties of the fluctuations. To 
do so, we will make use of the following standard results [23l |: 

Z = (2^) 3 V^tpI = (2tt) 3 [/3 2 - fi 2 .] 3 / 2 , (xiXj) = (A-% . (37) 

Then, it is easy to show that part of the energy is going into the fluctuations and that there is equipartition between 
the fluctuating parts of the magnetic energy and of the kinetic energy: 



(11 



One can also calculate the quantity of cross helicity going into the fluctuations: 



One should notice that there is no net magnetic helicity in the fluctuations because of the fact that A is strictly 
conserved. Then, the fractions of magnetic energy, cross helicity and kinetic energy going into the fluctuations are: 

^ - < U '' b '> - ^ M~\ (40) 



J B 2 dx JU-Bdx (3 2 - fi 2 c 
(u' 2 ) p 2 3/3 



M' 1 



where A4 = j B 2 efac is the magnetic energy of the coarsed-grained field. The first equation shows that there is an 
equal fraction of magnetic energy and cross helicity which goes in the fluctuations and the positivity of the magnetic 
energy requires: f3 2 > /j, 2 . Using this inequality and the second line, we can show that the fraction of kinetic energy 
going into the fluctuations is then bigger than that of the magnetic energy and cross helicity. This may gives some 
mathematical ground to the energy minimization procedure we used in section III CI 
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2. Rotating case 

The situation is changed when the other constant of motion are taken into account. We have: 

f b 

8H' m = 2 - dp dxdudb , (41) 
SL' = J ur Sp dxdudb . 



On the other hand. 



SL = J {8 AO + A SU) r 2 dx = J (U curP l 8B P + A SU) r 2 dx (42) 

(cui-r 1 ^ 2 ?/) ■ SB P + ASU r 2 ) dx 
(curl -1 (r 2 U) • bp + A u r 2 ) 8p dxdudb . 



Adding Lagrange multipliers —fjf m , —7 and —7' for H' m , L and L' respectively, we find that the expression (f 331) is 
multiplied by 

exp - i r u - 7 (curl -1 (r 2 U) ■ b P + A u r 2 )| . (43) 

The distribution of fluctuations is then still Gaussian and given by l|35() but now the mean-field equations are 

PV P + p c B P = , (44) 
PU + p c B + ~/r + jAr 2 = , 

f3B + 2p m A + p c U + ^ = , 
r 

[3B P + 2p m A P + p c XJp + cwr^rC'irA)] + 7 curV 1 (r 2 U) = . 

Taking the curl of the vectorial relations, we get the system (|13|l . 

Therefore, in the pesent case taking into account additional constant of motions, the relation between the velocity 
and the magnetic field is not linear anymore. The linearity is only valid for the poloidal component. The toroidal 
component obeys: 

f3(U+jr) = -fi c B~jAr 2 . (45) 

We can interprete U + 7' '/ 'fir as the relative velocity around a solid rotation £1 = — 7'//?. Indeed, 7' is the Lagrange 
multiplier for the angular momentum constraint. The non-trivial term responsible for the departure from linearity is 
— 7A?" 2 . Thus, the breaking of the proportionnality between the velocity and the magnetic field can be attributed to 
the conservation of the angular momentum in the Chandrasekhar model. This is an interesting feature because this 
conservation rule is likely to be more relevent in rapidly rotating objects. This may explain the dynamo saturation in 
rotating stars, where linearity between magnetic and velocity field is observed for slowly rotating stars and is broken 
for rotator faster than a certain limit (cf figure^). However, the non-proportionnality between velocity and magnetic 
field can also be due to additional conserved quantities such as those considered by Woltjer [l6| . 



IV. SUMMARY 

We have developped a statistical theory of axisymetric MHD equations generalizing the 2D approach by We 
derived the velocity and magnetic field distribution, and computed the corresponding equilibrium profiles for the 
mean flow. Like in the 2D case, the fluctuations around the mean field are found Gaussian, an universal feature 
connected to the conservation of the Casimirs under the coarse-graining. The equilibrium profiles are characterized 
by an alignment of the velocity and magnetic field, which is broken when the angular momentum conservation is taken 
into account. The statistical equilibrium profiles are found to correspond to profiles obtained under minimization of 
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FIG. 1: Magnetic field of stars (of late-type dwarfs) calculated from their X-ray emission, versus their rotation velocity. 



energy subject to the constraints. Thus, in the MHD case, in the presence of a coarse-graining (or a small viscosity), 
the energy is dissipated while the Helicity, the angular momentum and the Casimirs are approximately conserved 
(hydromagnetic selective decay). In particular, E — | J U 2 + B 2 dx ^ \ J (U 2 + B 2 )dx because part of energy goes 
into fine grained fluctuations Efi uct = E — E m j.. Therefore, the metaequilibrium state minimizes E at fixed /, H m , 
H c and L. This can be justified in the "classical case" (section IlII C 111 where we showed that the fraction of kinetic 
energy going into the fluctuating part of the fields was higher than that of the other quantities, namely the magnetic 
energy and the cross-helicity. The "rotating case" (section IIII C 21 requires more algebra and is left for further study. 

In contrast, in the 2D hydrodynamical case, the Casimirs are fragile quantities (because they are expressed as 
function of the vorticity which is not an integral quantity as the magnetic potential is) and thus are altered by 
the coarse-graining procedure. This is true in particular for a special class of Casimirs H = — j C(u>)dx., called 
ii-functions, constructed with a convex function C such that C" > 0. This leads to two very different behaviors of 
hydrodynamical turbulence compared to the hydromagnetic one. First, the H-functions calculated with the coarse- 
grained vorticity lJ increase with time while the circulation and energy are approximately conserved (hydrodynamic 
selective decay). Thus, the metaequilibrium state maximizes one of the -ff-functions at fixed E and T. For example, 
Chavanis and Sommeria 10] showed that in the limit of strong mixing (or for gaussian fluctuations), the quantity to 
maximize is minus the enstrophy, giving some mathematical basis to an (inviscid) "minimum enstrophy principle" . In 
this context, T2 = J uj 2 dx ^ J ui 2 dx because part of enstrophy goes into fine-grained fluctuations T fi uc t = T2 — T™ ■ 
However, for more general situations, the TJ-function that is maximized at metaequilibrium is non-universal and can 
take a wide diversity of forms as discussed by Chavanis 0] . Due to their resemblance with entropy functionals (they 
increase with time, one is maximum at metaequilibrium,...), and because they generally differ from the Boltzmann 
entropy Sb = — J u> In udx, the iJ-functions are sometimes called "generalized entropies" 0] . From the statistical 
mechanics point of view, there is an infinite number of constraints (depending on the micro scale fields) to take into 
acccount when deriving the Gibbs state. Consequently, the shape of the fluctuations is not universal. This is why the 
i/-function that is maximized at metaequilibrium is also non-universal. However, if the distribution of fluctuations is 
imposed by some external mechanism (e.g., a small-scale forcing) as suggested by Ellis et al. the functional S[uj] 
is now a well-determined functional determined by the Gibbs state and the prior vorticity distribution 0, 0] . 

Our computation can provide interesting insight regarding dynamo saturation. It is however limited by its neglect 
of dissipation and forcing mechanism. It would therefore be interesting to generalize this kind of approach to more 
realistic systems. In that case, the entropy might not be the relevent quantity anymore, but rather the turbulent 
transport, or the entropy production |2fij |. 
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APPENDIX: CURL OPERATORS 

Following Jordan and Turkington, we define 

curlB = (V x B) - eg (A.l) 
CurL4 = V x (A e e ) 

for any vector B and scalar A. It is straightforward to show that we have the following relations: 

curlCurl(-) = -rA*A (A.2) 



A curlB dx = I CurL4 • B dx (A.3) 
Setting A = Curl - B' and curlB = A' in the last identity, we get 

J Curl^B' A' dx = J B' • curh 1 / dx (A.4) 
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